How to count kinks: Prom the continuum to the lattice and back 
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We investigate the matching between (l+l)-dimensional nonlinear field theories coupled to an ex- 
ternal stochastic environment and their lattice simulations. In particular, we focus on how to obtain 
numerical results which are lattice-spacing independent, and on how to extract the correct effective 
potential which emerges from the simulations. As an application, we study the thermal production 
of kink-antikink pairs, obtaining a number density of pairs which is lattice-spacing independent and 
the effective barrier for pair production, i.e., the effective kink mass. 
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On 

^ , The possibility that the Universe underwent a series of symmetry breaking phase transitions during the earliest 
^ ' stages of its evolution has triggered a great deal of interest in the application of nonequilibrium statistical mechanics 
[t I , to cosmology. Of particular interest is the potential role that coherent field configurations, which arise from the 
interplay between nonlinearities and out-of-equilibrium conditions, could have played in shaping the earlier evolution 
and present-day structure of the Universe. Examples range from the nucleation of bubbles in the context of inflation 
CSJ ' and electroweak baryogenesis to the formation of topological defects |^ . 

^ , Given the relevance of the topic, and the obvious difficulties in performing experiments in a cosmological context, 
attempts to investigate the emergence of coherent field structures rely heavily on numerical simulations and possible 
analogies with condensed matter experiments 1^. Here we would like to focus on the former, namely, on numerical 
■ simulations designed to investigate the emergence of coherent structures in thermal field theories. An obvious limita- 
tion of such an approach is that, although field theories are continuous and usually formulated in an infinite volume, 
lattice simulations are discrete and finite, imposing both a maximum ("size of the box" L) and a minimum (lattice 
spacing Sx) wavelength that can be probed by the simulation. When the system is coupled to an external thermal (or 
quantum) bath, fiuctuations will be constrained within the allowed window of wavelengths, leading to discrepancies 
between the continuum formulation of the theory and its lattice simulations; the results will be dependent on the 
O • choice of lattice spacing. 
1^ , Parisi suggested that if proper counterterms were used, this depedence on lattice spacing could be attenuated 

This technique was implemented in a study of 2-dimensional nucleation by Alford and Gleiser Q. However, these 
^ studies still left open the question of how to match the lattice results to the correct continuum field theory. This is a 
crucial step if we want to test numerically certain predictions from field theories of relevance not only for cosmology but 
^ , also for condensed matter physics, such as decay rates of metastable states and the production of topological defects. 
Recently, Borrill and Gleiser (BG) have examined this question within the context of 2-dimensional critical phenomena 
. They have computed the counterterms needed to render the simulations indepedent of lattice spacing and have 
obtained a match between the simulations and the continuum field theory, valid within the one-loop approximation 
used in their approach. Inspired by their results, we decided to investigate the validity of this method within the 
context of topological defects. The results presented here should be relevant to numerical studies of the formation 
of topological defects and their comparison with experiments as well as to elucidating the general nature of 
the effective potential which emerges from coupling nonlinear field theories to a stochastic thermal (or quantum) 
background. 
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I. THE METHOD 



The Hamiltonian for a classical scalar field with potential Vo(0) and an environmental temperature T is, (with 
fcfi = c = 1) 
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Even though 1-dimensional field theories are free of ultra-violet divergences, the ultra-violet cutoff imposed by the 
lattice spacing will generate a finite contribution to the effective potential which must be taken into account if we 
are to obtain a proper match between the theory formulated in Eq. (|]) and its numerical simulation on a discrete 
lattice. If neglected, this contribution will compromise the measurement of physical quantities such as the density of 
kink-antikink pairs or the effective kink mass. However, before investigating the particular example of kink-antikink 
production, we present the method in its most general form. 

For classical, 1-dimensional finite-temperature field theories, the one-loop corrected effective potential is given by 
the momentum integral jsj 
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As mentioned before, the lattice spacing dx and the lattice size L introduce long and short momentum cutoffs 
A — tt/Sx and fc,„iii = 2TT/L, respectively. Lattice simulations are characterized by one dimensionless parameter, the 
number of degrees of freedom N = L/6x. For sufficiently large L one can neglect the effect of fcmin and integrate from 
to A. For Vq <C A^ (satisfied for sufficiently large A), the result can be expanded into 
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As is to be expected for a 1-dimensional system, the limit A — > oo exists and is well-behaved; there is no need 
for renormalization of ultra-violet divergences. However, the effective one-loop potential is lattice-spacing dependent 
through the explicit appearance of A, and so are the corresponding numerical simulations. In order to remove this 
dependence on 5x, we follow the renormalization procedure given by BG [||; it is irrelevant if the A-dependent terms 
are ultra-violet finite (d = 1) or infinite {d>2). In the lattice formulation of the theory, we add a (finite) counterterm 
to the tree-level potential Vq to remove the lattice-spacing dependence of the results. 
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There is an additional, A-independent, counterterm which was set to zero by an appropriate choice of renormalization 
scale. The lattice simulation then uses the corrected potential 



Matt(0) = 1^0 (0) 
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As we will show later in the context of kink-antikink pair production, this lattice formulation simulates the continuum 
limit to one loop as given by Eq. (|^). 

Note that the above treatment yields two novel results. First, that the use of Vlatt instead of Vq gets rid of the 
dependence of simulations on lattice spacing. [Of course, as 5x 0, Matt Vq. However, this limit is often not 
computationally efficient.] Previous works have explored the influence of a counterterm quadratic on lattice 
spacing. However, we note that for small enough 5x, the limit of interest here, our linear correction is dominant. 
Second, that the effective interactions that are simulated must be compared to the one-loop corrected potential 
Ml (0) of Eq. ^ ; once the lattice formulation is made independent of lattice spacing by the addition of the proper 
counterterm(s), it simulates, within its domain of validity, the thermally corrected one-loop effective potential. 



II. APPLICATION: THERMAL NUCLEATION OF KINK-ANTIKINK PAIRS 

As an application of the method discussed above we consider the symmetric double- well potential Vo((/)) — 
■| (0^ — (/ig) . The excitations of the associated quantum theory have a mass m = fiuj — ?i-\/2A0o- Thus, in order for 
the system to remain in the classical regime, the condition T ^ ?i-\/2A(/)o must hold. This constrains the dimensionless 
temperature Q = T/(-\/A(/'o) to be larger than \/2?i/(/)§. For <C Mk = ^8/9, where Mk is the dimensionless kink 
mass corresponding to the tree-level potential Vb |g] , we can expect to have only a dilute gas of kink-antikinks at 
thermal equilibrium. With these two conditions jointly satisfied, the system will also obey Mk = \/A(/)QMfc ^ to, 
indicating weak coupling. 

The corrected lattice potential is 
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simulations using Viatt will, in principle, match the continuum theory 

which has (shifted) minima at ±(/)iiii,i(r), with 4>uiin{T) < 0o- 

For the numerical simulations we introduce the dimensionless variables t = ^/Xipot, x — ^/X(j)QX, and (j) — 4'/4'Q- 
To keep the notation simple we will subsequently suppress the tilde. The field is prepared as (l){t — 0) = —1, and 
evolved in time according to a Langevin equation with white noise that incorporates the environmental temperature 
T through the fluctuation-dissipation theorem. The details of this and of the numerical implementation are laid out in 
The viscosity coefficient ij has been set to unity throughout this study. The time step is St — 0.05, and L = 2100. 
The heat bath takes a time Ai « 3 to achieve equipartition so that the energy per degree of freedom is E/N = T/2. 

III. RESULTS 

Ensemble average of field. For sufficiently low temperatures the simulated field will remain in the vicinity of 
the minimum (f> — — 0min(2^) for a very long time (compared to typical fiuctuation time-scales), until large-amplitude 
fluctuations drive portions of the space over the barrier at = and beyond. The subsequent evolution is then the 
formation of the first kink-antikink pair. True thermal equilibrium consists of reaching the final equilibrium kink- 
antikink density together with zero mean field. In a lose sense, this situation corresponds to symmetry restoration, 
although in one spatial dimension "symmetry restoration" will occur for any nonzero temperature; it is all a matter 
of time. 

As a first test of our procedure, we investigate the mean field value ^(i) = (1/i) J (j>{x,t)dx before the nuclcation 
of a kink-antikink pair, i.e., while the field is still well localized in the bottom of the well. In Fig. [l] we show the 
ensemble average of (j) (after 100 experiments) for different values of 6x, ranging from 1 down to 0.1, at T = 0.1. The 
simulations leading to the left graphs use the "bare" potential Vq, whereas the right graphs are produced employing 
Matt (Eq- Apart from a discrepancy for very coarse grids {6x = 1), where the resolution nears the correlation 
length, the average field value is clearly lattice-spacing independent when using Viatt, in contrast to the use of Vq. 

As discussed before, the average mean field value should correspond to the minimum —(t>min(T) of the effective 
potential. However, since we are only using a one-loop approximation, this agreement will get progressively worse 
as the temperature increases. For example, for T = 0.2, the discrepancy between the theoretical value, — (/)„iin(0.2), 
and the numerical result is 10%. For higher temperatures, we should not trust the one-loop approximation; other 
nonperturbative effects, such as subcritical fluctuations, too small in width and amplitude to emerge as a kink- 
antikink pair but still large enough to bring the average value of the field away from its one-loop value, will become 
important Thus, we restrict our investigation to temperatures safely within the limits of validity of the one-loop 
approximation. In a subsequent study, we intend to investigate the role of these nonperturbative effects. 

Density of kink-antikink pairs. Perhaps the most difficult task when counting the number of kink-antikink pairs 
that emerge during a simulation is the identification of what precisely is a kink-antikink pair at different temperatures. 
Typically, we can identify three "types" of fluctuations: i) small amplitude, perturbative fluctuations about one of the 
two minima of the potential; ii) full-blown kink-antikink pairs interpolating between the two minima of the potential; 
iii) nonperturbative fluctuations which have large amplitude but not quite large enough to satisfy the boundary 
conditions required for a kink-antikink pair. These latter fluctuations are usually dealt with by a smearing of the fleld 
over a certain length scale. Basically, one chooses a given smearing length AL which will be large enough to "iron 
out" these "undesirable" fluctuations but not too large that actual kink-antikink pairs are also ironed-out. In this 
study, a similar smoothing was implemented as a four-pole Butterworth low-pass filter of the field with a filter cutoff 
length AL. The filter removes fiuctuations with wavelengths smaller than AL. The choice of AL is, in a sense, more 
an art than a science, given our ignorance of how to handle these nonperturbative fluctuations. 

In Table 1 we show the number of pairs for different choices of fllter cutoff length and for different temperatures. We 
counted pairs by identifying the zeros of the flltered fleld. From Table 1 it is clear that as the temperature increases, 
the discrepancies in the count of pairs also increase. For this reason we only trust our data for fairly low temperatures. 
The problem is agravated by the fact that the "size" of the kink-antikink pair, i.e., the minimal separation between 
the two, not only changes due to dynamical effects, but also changes with temperature. Thus, choosing the fllter cutoff 
length to be too large may actually undercount the number of pairs. Choosing it too low may include nonpertubative 
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fluctuations as pairs. We chose AL = 3 in the present work, as this is the smallest "size" for a kink-antikink pair. 
In contrast, in the works by Alexander et al. a different method was adopted, that looked for zero-crossings for eight 
lattice units (they used Sx = 0.5) to the left and right of a zero crossing We have checked that our simulations 
reproduce the results of Alexander et al. if we: i) use the bare potential in the lattice simulations and ii) use a large 
filter cutoo length AL. Specifically, the number of pairs found with the bare potential for T = 0.2, 5x — 0.5 are: 
Up — 36, 30, and 27, for AL — 3, 5, and 7, respectively. Alexander et al. found (for our lattice length) Up = 25. 
Comparing these with Table 1, it is clear that the differences between our results and those of Alexander et al. come 
from using a different potential in the simulations, viz. a corrected vs. an uncorrected potential. 

We believe that at this point it is fair to say that the "smearing issue" remains unresolved, at least for temperatures 
T > 0.25 or so. We intend to address the issue of how to deal with these nonperturbative effects in a forthcoming 
publication. In any case, the focus of the present work is mostly on how to achieve a lattice-independent count, 
irrespective of the particular method used for identifying the kink-antikink pairs. 

Fig. 1^ compares measurements of the kink-antikink pair density (half the number of zeros of the filtered field), 
ensemble-averaged over 100 experiments, for different lattice spacings. Again it is clear from the graphs on the left 
that using the tree-level potential Vq in the simulations causes the results to be dependent on 6x, whereas the addition 
of the finite counterterm removes this problem quite efficiently; both diagrams of Fig. ^ contain four graphs each, 
although the graphs on the right are almost indistinguishable. Unless the properly corrected potential is used in the 
lattice simulations, the measured number density of topological defects is sensitive to the lattice spacing. One must 
be careful when counting kinks, especially for large lattice spacings, say Sx = 0.25 or larger. 

The next step is to extract the correct continuum theory from the lattice simulations. What theory is the lattice 
simulating? Most previous simulations of thermal nucleation of kink-antikink pairs have overlooked this problem. 
Although a temperature-dependent kink mass was conjectured in the works of reference ]TT[ |, not much has been done 
to understand its origin or its value. One way of addressing it is by comparing the numerically measured kink mass 
with its theoretical prediction. It has been found that the measured mass was smaller than the theoretical prediction 
by a factor ranging from 25% to 45% ||ll|,|l^, a disturbing result. This has been attributed to several effects, such as 
the finite size of the lattice, the finite size of the kinks, and phonon dressing effects due to the lattice discretization 
p3[ . We will show that this problem is rooted in the incorrect matching between theory and numerical simulations. 
In the works by Alexander et al. a beautiful agreement between the low temperature limit and a T = WKB 
approximation was obtained, as well as between high temperatures and a double Gaussian nonperturbative method 
pO| . Our method is effective precisely between these two regimes, and could be interpreted as a T-dependent WKB 
approximation obtained naturally from the inclusion of counterterms. 

One should expect the equilibrium kink-antikink pair density to follow the proportionality 
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where Mk is the kink mass, given by 
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and (/)/c(x) is the kink solution to the equation of motion. Note that we left the potential V{(l)k) unspecified. If we 
use the tree-level potential, Vb((/)fc), we obtain the well-known result Mk = \/8A/9(/)q. Or, in dimcnsionless variables, 
Mk = •\/8/9. One can extract the numerical value of Mk by measuring the pair density and plotting the results in a 
logarithmic scale, as in Ref. |l^. The result should be a straight line with negative slope —Mk = — \/8/9. However, 
as mentioned above, the measured slope was found to be about — 0.70A/fc. The reason for the discrepancy is that 
the potential which should be used when comparing theory and simulation is not the tree- level potential Vo{4>) but 
the effective potential Vil (</>). Thus, one must compute the effective kink mass Mk{T) using the corrected potential 
yn^{4') s-nd then compare the results with the numerical simulations. 



The effective kink mass can be found using the equation of motion and the real part of Vn^ |15 
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This integration can easily be carried out numerically. In Fig. |^ we plot the ratio Mk{T) / Mk vs. the dimcnsionless 
temperature, 8. Of course, for T = 0, Mk{0) / Mk = 1. As the temperature increases, the effective kink mass decreases. 
The points represent the kink mass extracted from the numerical simulations, while the error bars were obtained by 
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propagating the standard deviation of the ensemble average. It is quite clear that the effective kink mass tracks the 
numerical values quite well. In fact, within the validity of our approximations, the "averaged" value for the effective 
kink mass is 0.75Mk. Also, since the mass extracted from the simulations depends on the filter cutoff length AL, the 
reasonable agreement between theory and numerical experiment offers indirect support for our choice of AL = 3. For 
very small and very large temperatures the theory fails to track the numerical data. At large temperatures > 0.25, 
the one-loop approximation breaks down, while for low temperatures & < 0.12, the large pair nucleation time-scale 
precludes a proper statistical analysis (not enough experiments). However, the conclusion is quite clear: by controlling 
the dependence on lattice spacing of the simulations we were able, within the validity of our approximations, to obtain 
the correct effective potential that should be used when comparing theory and numerical experiment. 
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FIG. 1. Average field value for T = 0.1 using the tree-level potential, left, and the corrected potential, right. 




FIG. 2. Density of kink-antikinks (half of density of zeros), for T = 0.2 and 6x = 1, 0.5, 0.2, and 0.1. 
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FIG. 3. The ratio of the effective kink mass, Mfc(T), to the uncorrected kink mass, Mk, vs. the temperature. 



AL 


T = 0.15 


T = 


= 0.20 


T = 0.25 


3 


10.3 lb 0.2 


39.0 


lb 0.5 


75.6 lb 0.6 


5 


8.9 ± 0.2 


32.7 


lb 0.5 


62.0 ± 0.5 


7 


8.4 ±0.2 


29.8 


lb 0.4 


54.6 lb 0.5 



Table 1: Number of kink-antikink pairs for different choices of the filter cutoff length AL and T. 
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